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Abstract 

We discuss the applicability of the perturbation theory in electrodynamic prob- 
lems where the local Leontovich (the impedance) boundary conditions are used to 
calculate the ohmic losses at the metallic surface. As an example, we examine a 
periodic grating formed from semi- infinite rectangular plates exposed to the s- 
polarized electromagnetic wave. Two different ways for calculation of the ohmic 
losses are presented: (i) the calculation of the reflection coefficient obtained with 
the aid of the perturbation theory (the impedance is the small parameter) and (ii) 
the direct calculation of the energy flux through the metallic surface, when to get 
the answer only the tangential magnetic field at the surface of a perfect conductor 
of the same geometry has to be known. The results (i) and (ii) differ noticeably. 
The same difficulty is inherent in all the problems where the metallic surface has 
rectangular grooves. We show that the standard first order perturbation theory is 
not applicable since beginning from a number n even the first corrections to the 
modal functions <p n used to calculate the fields, are of the same order as the zero 
order modal functions (the impedance is equal to zero). Basing on the energy con- 
servation law we show that the accurate value for the ohmic losses is obtained with 
the aid of the approach (ii) . 



1 INTRODUCTION 



I _ 

*t3 situation arising when the perturbation theory is used for calculation of the ohmic losses 
G at the highly conducting metallic surfaces with sharp edges. As an example we examine 
(j the one-dimension metallic periodic grating formed from semi-infinite rectangular plates 
j>! (see Fig.l) that is illuminated by s-polarized light. 

All the calculations are performed in the framework of the local impedance boundary 



conditions (the Leontovich boundary conditions) applicability |T], |2|]. In other words, we 
assume that 

<5<6<A, (1) 

where 5 is the penetration depth of an electromagnetic wave into the metal, b is the 
characteristic size of the surface contour (in our case, the period of the grating) and 
A = 27rc/uj is the vacuum wavelength, u> is the frequency of the incident wave. 

This allows us to restrict ourselves to solution of the Maxwell equations in the region 
external with respect to a metal. At the metallic surface the boundary conditions for the 
fields are 

E t = C[n,H t ], (2) 
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E t and H t are tangential components of electric and magnetic vectors at the metallic 
surface, n is the unit external normal vector to the surface. The two-dimensional tensor ( 
is the surface impedance tensor. The real and imaginary parts of the impedance define the 
dissipated energy and the phase shift of the reflected electromagnetic wave respectively 
(see, for example, In what follows we suppose that in Eqs.(2) the impedance ( is 
an ordinary multiplying operator, and it is the same as the impedance of a perfectly flat 
surface of the same metal. 

Of course, the boundary conditions (2) are not rigorous ones. In general the relation 
between E 4 and H t is nonlocal. Moreover, when the metallic surface is not flat, even 
the local terms in the expression for the components of the surface impedance tensor 
have corrections depending on the shape of the metallic surface (see, for example, || £§.) 
However, when the inequalities (1) are fulfilled, the leading term in the expression for 
( is equal to the impedance of the perfectly flat metallic surface. Having in mind to 
calculate the leading term in the expression for the ohmic losses we restrict ourselves to 
these approximation only. 

In what follows we assume that the metal is an isotropic one. Then the surface 
impedance tensor is ( Q/ 3 = C<W- Under the conditions of normal skin effect C = 
^Jujjsna(l — i), where a is the conductivity of the metal. Under the conditions of anoma- 
lous skin effect the expression for the surface impedance is more complex (see, for example, 
H). In our calculations we do not specify the character of skin effect. We only take into 
account that in the order of magnitude |£| ~ 5/X. For good metals this ratio is extremely 
small. 

Note, the local impedance boundary conditions (2) are valid up to the terms of the 
order of |£| only. Consequently, when calculating the fields the terms of the higher order 
(for example, of the order 6|C|/A, |C| 2 and so on) must be omitted as far as they exceed 
the accuracy allowed by the boundary conditions. 

Strictly speaking, inequalities (1) provide the possibility to use the local boundary 
conditions, Eqs.(2), when calculating the fields above the metal only if the surface of the 
metal is rather smooth. The presence of sharp edges makes this possibility questionable. 
Of course, in real life all the sharp edges are rounded ones. Undoubtedly, the boundary 
conditions (2) can be used only if the smoothing characteristic size is greater than 5 too. 
However, it is rather difficult to take the smoothing regions into account accurately even 
for very simple surface profiles. Therefore rather often Eqs.(2) are used, assuming that a 
smoothing changes the results only slightly. 

It must be noted that the grating depicted in Fig.l is infinitely long in the direction. 
At first it seems that the boundary conditions (2) are not valid in this case. To justify 
ourselves we take into account that in the s-polarization state the fields decay exponen- 
tially inside the grooves when the distance from the throat of the groove increases. As 
a result the fields in the vicinity of infinitely deep grooves and the grooves of finite, but 
rather large depth h (h > ma and m is of the order of some units) are practically the 
same. For such grooves the smoothing of the contour allowing us to use the boundary 
conditions (2) can be done rather easily. 

After all these preliminary remarks we are faced with a rigorously set mathematical 
problem: suppose that a plane wave E l = E e~^ fca:3+a; ^ (k = uj/c) is normally incident 
upon the surface depicted in Fig.l. At the metallic surface we have the boundary con- 
ditions (2) that are valid up to the terms of the order \C\. We impose an outgoing wave 
or exponentially decaying boundary conditions on the reflected fields. Our problem is to 
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define the ohmic losses. In other words, we need to calculate the energy flux through the 
surface of the metal. 

Let us introduce the time-averaged Pointing vector in the region outside the metal 

S = (c/8vr)Re[EH*] (3) 

that is the density of the time-averaged energy flux. According to the energy conservation 
law 

divS = 0. (4) 
The consequence of the last equation is the equality of the integrals 

/ Sdf = / (Sn)d/, (5) 

where the integral in the left-hand side is taken over the surface infinitely distant from 
the metal. The integral in the right-hand side is taken over the surface of the metal. It 
defines the energy flux penetrating through the metallic surface. 

It is well known that infinitely far from the metallic surface (x^ ^> A) the only nonzero 
component of the Pointing vector is expressed in terms of the reflection coefficient R: 
S3 = — (c/87r)|Eo| 2 (l —R). Consequently, when the reflection coefficient is known, we can 
calculate the ohmic losses with the aid of the left-hand side of Eq.(5). 

On the other hand we can calculate the energy flux through the surface of the metal 
directly. With regard to the boundary conditions (2) at the metallic surface the normal 
component of the time-averaged Pointing vector is S n = ^rRe£|H t | 2 . The last expression 
can be simplified. 

Since \(\ <C 1, the magnetic vector H t at the surface of a good metal is approximately 
the same as the magnetic vector at the surface of the perfect conductor (the conduc- 
tivity is infinitely large and, consequently, the impedance ( = 0) of the same geometry 

0. Then within the allowed accuracy we must rewrite S n as S n = ^ReC|Hj°' ) | . 

Let us introduce the energy Qd dissipating in the metal per unit length in the 22 
direction and per one period. According to the previous argumentation and making use 
of Eq. (5) we can write two equivalent expressions for Qd 

Q d = ^ReC / |H| 0) | 2 d/; (6.0) 
Sir Jl 

the integral is taken over the contour of the surface relating to one period, and 

Q d = 2b\E \ 2 ^[l-R\. (6.6) 

Note, when Eq.(6a) is used, to calculate Qd we need to know the magnetic vector 
at the surface of the perfect conductor only. On the other hand, when we make use of 
Eq.(6.b) we need to know the reflection coefficient R. To calculate R we need the solution 
of the complete electrodynamic problem in the region over the real metal. 

A lot of authors calculated the fields above perfectly conducting gratings (see, for 
example, f6j). It can be done with the aid of different methods @. For example, the fields 
inside the grooves are presented as series of modal functions that are the solutions of the 
Maxwell equations. The representation incorporates the boundary conditions for perfect 
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conductors. Above the grooves the usual Rayleigh plane-wave representation is adopted. 
As a result a matrix equation for the modal amplitudes is obtained. This method makes 
use of the standard Fourier analysis. There are no doubts with respect to the results 
obtained with the aid of this method. The other way is based on the solution of the 
Wiener- Hopf equation. Both methods give the same result. 

Applying the local impedance boundary conditions, we can use the same methods 
when calculating the fields above highly conducting metallic gratings. However, to be 
sure that we do not fall outside the limits of the boundary conditions (2) applicability, 
when calculating we have to exploit the perturbation theory based on the inequalities 
(1). In what follows we show that in this problem the standard perturbation theory is 
not applicable. More accurately, we show that in our problem the dissipated energy Qd 
calculated with the aid of Eq.(6.a) differs from the result obtained with the aid of Eq.(6.b) 
after calculation of the reflection coefficient R in the framework of the perturbation theory. 
This discrepancy leads to violation of the energy conservation low. 

We would like to stress once more, the field entering Eq.(6.a) can be calculated 
with high accuracy. Then, in the framework of our approximation Eq.(6.a) provides the 
correct result for the ohmic losses. Next, as it is shown below, the method of calculation 
of the fields above the real metal (the impedance ( 7^ 0) is much alike the one used when 
calculating the fields above the perfect grating. However, to calculate the difference 1 — R 
up to the terms of the order of (, we have to use the perturbation theory with respect to 
the small impedance (. Consequently, the source of the obtained discrepancy can be the 
perturbation theory only. 

We faced this problem when calculating the ohmic losses on the lamellar grating with 
finite depth of the grooves exposed to the s-polarized light ||. We found out that for 
rather deep grooves our result was approximately twenty percents less than the result 
obtained by L.A.Vainshtein, S.M.Zhurav and A.I.Sukov in |J. In the same frequency 
region defined by Eqs.(l) in the framework of the impedance boundary conditions (2) 
these authors examined the one-dimensional grating with infinitely deep grooves and very 
thin plates ((& — a)/b <^ 1). They based their solution on the application of the Wiener- 
Hopf equation and calculated the reflection coefficient. Then taking into account that 
according to Eqs.(l) the parameter \(\/kb <C 1, they used the perturbation theory and 
calculated the difference 1 — R entering Eq.(6.b) up to the terms of the order of (. 

In the s-polarization state there are no physical reasons of different results for very 
deep and infinitely deep grooves. Therefore, trying out to verify our approach, for the 
grating with infinitely deep grooves we reproduced the perturbation theory calculation 
of with the aid of the modal function method adopted in our calculations in [JSJ. The 
results obtained with the aid of the perturbation theory were the same as in the work of 
Vainshtein and his co-authors. Then we repeated the calculation with the aid of Eq.(6.a). 
Comparing the results we found the same twenty percents difference (see Fig. 2). 

The result of this analysis is presented in this publication. From our point of view 
our calculations show that in this problem the application of the standard perturbation 
theory provide a regular error when calculating the fields above the grating. In what 
follows we show the place where the standard perturbation theory failed. 

It was not the goal of this work to develop a regular perturbation theory allowing us 
to calculate the fields above lamellar conducting gratings correctly. The main results of 
our investigation can be formulated as follows: 

i. For highly conducting lamellar gratings the correct result for ohmic losses Qd is given 
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by Eq.(6.a). 

ii. The coincidence of the result for Qd obtained with the aid of Eq.(6.a) and Eq.(6.b) 
can be used as a proof of correct calculation of the fields above the metallic surface. 

The organization of this paper is as follows. In Section 2 we present the modal function 
approach allowing us to calculate the fields above the metallic grating. In Section 3 
we describe the results obtained in the framework of the standard perturbation theory. 
In Section 4 we calculate the ohmic losses and discuss the reasons why the standard 
perturbation theory fails. Concluding remarks are given in Section 5. In APPENDIX we 
show that for an arbitrary value of kb < 1, the formulae used when calculating the fields 
above the infinitely conducting grating provide the total reflection of the incident wave 
(the reflection coefficient R = 1). 



2 CALCULATION OF THE FIELDS 

Let the one-dimensional periodic metallic grating depicted in Fig.l be exposed to the 
s-polarized normally incident electromagnetic wave E J = (0, E , 0)exp[— ikx 3 \. For sim- 
plicity we assume that on the flat parts xs = of the surface the impedance is equal to 
zero. On the inner surfaces of the semi-infinite grooves the boundary conditions (2) are 
fulfilled. 

The inequalities (1) provide us with three small parameters. Namely, 

|C| ~ T < !> kb < 1 and T ~ T7 < !■ ( 7 ) 
A o kb 

No other restrictions are imposed on the ratio a/b < 1. 

At first, we will not pay attention to inequalities (7). In what follows we obtain the 
solution of the problem without using the perturbation theory. We assume only that 
kb < 1. 

Let the superscripts + and — denote the fields in the half-space x 3 > and in the 
central groove respectively. With regard to the surface periodicity we adopt the represen- 
tation 

E+(x u x 3 ) = E | e~ lkX3 + jr e +e 4(m?Xl/fc+a ^ 3) 1 , (8) 



for the electric field above the metal. In Eq.(8) ao = k and a q = iy(irq/b) 2 — k 2 , if 
q 7^ 0. Since on the plain parts x 3 = of the surface the impedance is equal to zero, 
i?2~(a;i,0) = 0, when a < \x\\ <b. 

Inside the central groove (\%i\ < a, x 3 < 0) we present E~(xi, x 3 ) as a series of modal 
functions 4> n (xi,x 3 ) that are the solutions of the Maxwell equations incorporating the 
boundary conditions (2) on the vertical sides of the groove?]. Let B n be the coefficients of 
this series. 

Separating the variables we write the modal functions (p n (xi,xs) as 



4> n (x 1 , X 3 ) = ^(xaVnO&l); ^n{x 3 ) = exp 



- (ka) 2 



(9.1) 



^^An analogous method was repeatedly used when calculating the fields above infinitely conducting 
rectangular gratings (see, for example, M) 
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It is easy to see that the functions (p n (xi) are given by different expressions for even and 
odd numbers n: 

sin(7 2n xi/a) cos(7 2n _iXi/a) 

f2n\xi) = ; tpm-iixi) = : (9.2) 

COS72n sm7 2 „-i 

with 7 2n and j2n-i being the solutions of dispersion equations 

tan7 2 „ i( cot7 2n _i i( 



l2n ka 72n-i ka 



(9.3) 



The set of modal functions 0„(xi,x 3 ) has a lot of interesting properties. We do 
not discuss them here, but refer the reader to the review work of Yu.I.Lyubarskii ||10|| . 
Note, the functions <p n (xi) are not orthogonal ones since their scalar product (y? m ,y? n ) = 
(l/2a) /" a dxi(pm(xi)(p n (xi) (the star means the complex conjugate) does not equal to zero 
when m 7^ n. However, the functions (p n (xi) with different subscripts are orthogonal: 

1 f a 

— J dXiV2m(^l)v 2 n(Xl) = ImmSmn, ( 10 -1) 



where 



hm,2m — 7^7 \ o ^ T~ f ' hm-l,2m-\ — 7TT \ "T^ ^ T~ \ • (10.2) 

2b lcos / 7 2m ka \ 2b Ism 72m- 1 ka 



We used this property of the functions (p n (xi) when deriving the matrix equation for the 
modal coefficients (see below Eq.(13)). 

Now with regard to Eqs.(9) we write the field E~(xi,x 3 ) (\xi\ < a, x 3 < 0) as 

{00 00 ^ 

B 2 n1p2n(x3)<P2n(xi) + B 2n-1^2n-l(x3)<P2n-l(xi) > ■ (11) 
n=0 n=l ) 

Applying the boundary conditions on the plateaus x 3 = and the continuity conditions 
across the central slit, we obtain two linked matrix equations allowing us to determine 
both the amplitudes e+ and the modal amplitudes B n . In place of the amplitudes B n it 
is convenient to introduce the coefficients X n , 

i{kb)X n = ^B n . (12) 

7T 

We eliminate the amplitudes e^ + - ) from the set of matrix equations to obtain the 
matrix equation for the coefficients X n . There are two separate matrix equations for 
the coefficients X 2n and X 2n -i- The even coefficients X 2n satisfy a homogeneous matrix 
equation. Consequently, all the coefficients X 2n = 0. For the odd coefficients X 2n _i our 
calculations results in the following matrix equation: 



oo 



2 



; \b \ • it 

X 2m -1^2m-l,2l-l + [ 

m=l 

If we set fi = 2a/b, the matrix A2 m _i ) 2z-i m Eq.(13) is 



]T X2 m _iA 2m _i,2i-i + ( ^ ] ( 7T— ) V 1 " (ka/l2i-i) 2 l2i-i,2i-iX 2l ^ = 2W a -i. (13) 

2a I \ 2y 2 i-i J 



A 2m - ll2 z-i = J2 W 1 - {kb/mfD2m-i,2i-i{q) - i-(kb)W 2l ^W 2m ^; (14.1) 

9=1 71 
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C2«-i,g stands for the integral 

(^jr 1 ) Ca-i,, = Y b j^^ie-^'ViW- (14.3) 
The explicit form of C 2 i-i >q is 

//!7r\ cosw, - WgSin[^cot7 2 /_i/72«-i] vrg/x 

^-^-(tJ ^ - 72 Vi ' Uq -~T' { ] 

and hi-1,21-1 is defined by Eq.(lO). 

We can calculate the coefficients X 2n ~i solving numerically the infinite set of equations 
(13). To calculate the reflection coefficient R we take into account that since kb < 1, there 
is only one reflected wave. In accordance with Eq.(8) its dimensionless amplitude is 
and e+ with q ^ represent evanescent waves. Omitting intermediate calculations, we 
present the expression for eg in terms of the coefficients X 2 i-i- 



e + = _i -i^)^^!^. (15) 



Then i? = |4|. 

It seems, that in this way we solve the problem and can calculate the ohmic losses sub- 
stituting the found reflection coefficient into Eq.(6.b). However, as it was aforementioned, 
the local impedance boundary conditions (2) are true up to the terms of the order of ( 
only. Consequently, to obtain a reliable answer for the ohmic losses we need to extract 
the terms of the order of ( from the expression for 1 — R. And this is just the place where 
the perturbation theory has to be used. 

3 THE RESULT OF THE FIRST ORDER PER- 
TURBATION THEORY 

The zero term of our perturbation theory corresponds to £ = 0. We use the superscript 
(0) to denote that £ = 0. Then in accordance with Eqs.(9) we have 

7& ) -i = 7r(2W 2 ~ 1) ; ^S_i(x 1 ) = (-l) B cos(7r(2n-l)x 1 /2a). (16) 

The functions v4n-i( a; i) are cosines, so they can be used as basic functions for the standard 
Fourier analysis. 

Next, from Eq.(9.3) it follows that the eigenvalues 721-1 depend not on ( itself, but on 
the parameter e = (/(ha). Since ka <C 1, it is evident that \(\ <C \e\ <C 1. With regard 
to Eq.(9.3) we obtain that up to the terms of the order of e 

721-1 = 7^(1 -ie). (17) 

Moreover, the impedance ( enters the matrix equation (13) only through the dependence 
on 722-1. Thus, from Eq.(13) it follows that the coefficients X 2n -i depend not on all the 
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three small parameters listed in Eq.(7) separately, but on the two small parameters only. 
They are kb and e. When e = the coefficients X^_ 1 (kb) define the fields above the 
infinitely conducting grating. With regard to Eqs.(13) and Eqs.(14) it is easy to see that 
the matrix equation for X^-iikb) is 



E^-iAgL^ + Vi - (kb^(2i - lyf^ = (i8.i) 

where 

OO 1 

^2m-i,2Z-i(^) = E W 1 - (kb/vqfDfl^M - Akb)- (18.2) 

and 

= = to,)2 - (2m - 1)2 (18 ' 3) 

uj q is defined in Eq.(14.4). 

Suppose that we know the solution of Eq.(18.1). To take into account the finite 
conductivity of the metal, we need to make use of the perturbation theory when solving 
the matrix equation (13) and calculate the first corrections to X^-i(kb) in the small 
parameter e. 

Note, according to Eq.(ll) and the definition of the coefficients X 2n -i, Eq.(12), the 
same corrections provide the electric field E~ up to the terms of the order of (. Since 
the starting point of these calculations, the local impedance boundary conditions (2), 
do not allow us to calculate the fields up to the terms higher than (, when calculating 
the coefficients X 2n _i the terms of the order higher than e (for example, of the order 
of () must be omitted within the accuracy of the local impedance boundary conditions 
applicability. 

With the aid of Eq.(17) we expand all the coefficients in Eq.(13) up to the terms of 
the order of e. Then with regard to Eqs.(18) after tedious, but straightforward algebra 
we obtain that up to the terms of the order of e the solution of Eq.(13) has the form 

X 2n . 1 (kb; e) = Xiti(kb) + 2ie(X$>_ 1 (0) ~ V^_i), (19) 
where the coefficients V 2n -i satisfy the matrix equation 

l/7T\ 2 V^-l 
V2m-l^ 2m -l,2l-l{V) + ' 

m=l 



E V 2m .^L^M + - o 9HTT = S *£Uo)S 2m -i, 2 *-i; (20.1) 

71=1 A* VZ/ U L m=l 



the matrix A^-i 2j-i(0) * s defined by Eq.(18.3) when kb = 0, the matrix SWn-i^-i i 

00 1 /7T\ 2 1 

S 2m -l,2l-l = E(^) T 2m-l,2Z-l(<?) + -(^2 ) ^737^' ( 20 - 2 ) 

and 



is 



T 2m -i,2i-i{<i) 



d m {q)di(q) 



. ({2m -l) 2 (2/-l) 2 \ . o 

cos 2uj q j 1 + uj q sin 2u q 

\ d m di I 



(20.3) 

Again uj q = nqfi/2 and d m is defined in Eq.(18.3). 

If in consecutive order we solve Eqs.(18) and (20), we obtain the coefficients X 2n ^i(kb, e) 
up to the terms linear in the parameter e. 
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4 CALCULATION OF OHMIC LOSSES 

At first, let us calculate the dissipated energy Q d with the aid of Eq.(6.a). We mark the 
result of this calculation with the superscript (sur). We remind that we set £ = at the 
flat parts £3 = of our surface. Then in Eq.(6.a) the integration is performed over the 
vertical facets of the central groove, where the tangential magnetic vector is H 3 °\±a, x 3 ) 
(x 3 < 0). Now we can rewrite Eq.(6.a) as 

Jsur) _ C [0 (0) ,,2 



Qr> = ^ReC J ^ \H^(a,x 3 )\ dx 3 ; (21) 

We take into account that since the right-hand side of Eq.(21) is proportional to (, within 
the accuracy of our approach the field H^\a : x 3 ) in the integrand has to be known only 
up to the terms independent of the small parameter kb. Then with regard to the Maxwell 
equations and equations (11) and (12) (we remind that all the coefficients I?2n = 0) we 
have 

00 

Hi°\a,x 3 )\ kb=0 = — £oE X 2?li(0)exp[7r(2Z - l)x 3 /2a], (22) 
V J=t 



and, consequently, 



= cb\E \ 2 ^ £ (23) 



8// , =1 l + m-1 



Now let us calculate the energy Q d with the aid of Eq.(6.b). We mark this result with 
the superscript (00). Since the reflection coefficient R = | Cq" | 2 , to calculate we use 
Eq.(15). Then 



(oo ) _ cb\E Q , 
y d - An V, 

a 00 /u\ 2 00 

Q = i(kb)^J2( X 2l-l W *-l- X *2l-l W 2l-l) + (M) 2 - E 4-1^-1^,-1^-1. 
71 1=1 V7F/ l,m=l 

(24) 

When calculating Q d °°^ it is necessary to have in mind, that the dissipated energy 
has to be equal to zero for a perfectly conducting surface. In other words, when ( = 
and, consequently, X 2 i-\ = X^}_i(kb), Q d °°^ vanishes for an arbitrary value of kb < 1. 

In APPENDIX we show that this equality is fulfilled for X^iL^kb) that are solutions of 
Eq.(18.1). 

As a result, substituting the expression (19) in Eq.(24) we obtain Q d °°^ in terms of the 
coefficients X^l-^O) and Vy-i- 



2 



71 1=1 — J-J 

This equation is valid up to the terms of the order of £ 

Our next step was numerical calculation of the coefficients^ X^l^O). With this result 
in hand we evaluated the right-hand side of the matrix equation (20.1) and calculated 

2 The computational procedure used in the present work when calculating the coefficients X^l^O) is 
just the same as the one used in ||. So we refer the readers to the work ||. 
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the coefficients V21-1. Finally, we calculated numerically Qj T and with the aid of 
Eq.(23) and Eq.(25) respectively. 

In Fig. 2 we plot the dimensionless values % = 7tQ £° / cb\Eo\ 2 Re( and q s = ttQ / cb\Eo 
as functions of the ratio a/b. We see that for all a/b < 1 the results for and 
differ noticeably. 

In Fig. 2 we also display the results for the ohmic losses obtained by L.A.Vainshtein 
and his co-authors in ||. As it was mentioned in INTRODUCTION, they examined the 
same grating formed from thin plates ((b — a)/b <C 1). In || the way of calculation was 
much alike the one used in the present work when calculating Q^°\ however some other 
modal functions were used, and the problem was reduced to solution of the Wiener-Hopf 
equation. Supposing that the impedance of the flat parts x% = of the metallic surface 
was equal to zero, the authors of obtained an elegant analytic formulae for the leading 
term of the ohmic losses 

gf) = ^Re C ln(-^), (b-a)/b«l. (26) 
7i z \b — a) 

(The superscript (V) indicates that this is the result of ||.) 

We see, that our results for are in good agreement with the values ■ This 
confirms our calculations made in the framework of the perturbation theory. However, 
neither , nor coincide with that is the true energy flux penetrating through 

metallic surface. Consequently, we conclude that the obtained result for (the same 
as Q4) is improper, since it leads to violation of the energy conservation law. 



4.1 Inefficiency of the standard perturbation theory 

We think that the source of this discrepancy is application of the standard perturbation 
theory. The point is that if ( 7^ the set of modal functions (f>2n-l( x l, ^3) (see Eqs.(9)), 
used when calculating the reflection coefficient, being complete and minimal still does not 
form a stable basis. As a result expansions in this set can be very sensitive to perturbations 
(see, for example, [|T0[] ). 

To show this, let us note that according to Eqs.(16) the zero order eigenvalues 72°- 1 
are equidistant numbers: 72° + i — ifn-i = 71 ■ However, the distance between 72,1-1 defined 
by Eq.(17) and 72°_i increases linearly when n increases. If, for example, ( = — i) 
(that is true for homogeneous solids under the conditions of normal skin effect), we have 



1 (°) 1 

\l2n-l -72n-ll * 



P s/2 



(2n - 1), (3 = (o/ka. (27) 



Now let us calculate the eigenfunctions f2n-i( x i) substituting 72 n -i from Eq.(17) in 
Eq.(9.2). The result of this calculation is 



<£2n-i(zi) = (-1) 



n cosh/foff-i^ cos(l + (3hiTif - * sinh^l^ sin(l + (3)4°lif 
cos /?72° } _i cosh /372?-i - i sin /fy^i sinh fi^n-i 



(28) 

We see, that y?2n-i(^i) is a complex function, whose real and imaginary parts are of 
the same order. The difference between <^ 2 n-i(^i) arid y4n-i( x i) is small only if n <C 
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l/n/3. When n ^> l/n(3 the eigenfunctions (p2 n -i( x i) do not resemble the non-perturbed 
eigenfunctions y4n-i i x i) a t all- In Fig.3 for x\ = we plot the ratios Rey22n-i(0)/y4ri-i(0) 
and Im<^2n-i(0)/v32n-i(0) fo r two values of f3: (3 = 0.05 and (3 = 0.01. 

We must also take into account that in the s-polarization state the strength of the 
magnetic field in the immediate vicinity of a rectangular two-dimensional wedge increases 
significantly (see, for example, When the fields are represented by the expressions of 
the type (11), to describe the fields near the corner points correctly we need to provide 
a special asymptotic behavior of the coefficients X 2n _i when n — > 00 (see 0). For ( = 
the coefficients -X^n-i ~ 1/ (2n — l) 2 ^ 3 when n — > 00 (see ||) 

With regard to all the aforementioned arguments we take for granted that in the 
s-polarization state the perturbation theory presented in Section 3 is not applicable 
when the local impedance boundary conditions are used to calculate the fields above 
one-dimensional periodic rectangular metallic gratings. 



In this section we summarize the results of our analysis. We would like to remind that 
here we examined the case of the s-polarization only. We put aside the p-polarization 
state. In the end of this section we'll add some remarks relating to the p-polarization 
state. 

From our point of view, the main result of this publication is the demonstration that 
the standard perturbation theory can fail if in the framework of the local impedance 
boundary conditions one try to use it when calculating the electromagnetic fields above 
periodic metallic gratings. We examined the case of rectangular infinitely deep grooves. 
However, it is quite possible that the same problems arise when the shape of the grooves 
is, for example, a triangular one, or, more general, when the surface profile has sharp 
edges. 

Our calculations showed that for rough surfaces the results for the reflection coefficients 
obtained with the aid of the perturbation theory (the impedance ( is a small parameter) 
are not always reliable. On the other hand, there are no doubts that the dissipated energy 
calculated with the aid of Eq.(6.a) defines the ohmic losses accurately. Thus, we can define 
the reflection coefficient comparing equations (6. a) and (6.b). Then 



where the tangential magnetic field H t is calculated up to the terms independent of kb. 
Let us remind that our results were obtained in the frequency region consistent with the 
inequalities (1). 

On the other hand, this approach allows us to calculate the reflection coefficient only. 
The fields in the close vicinity of a rough metallic surface (at the distances d that are of 
the order of the period of the surface structure where the evanescent waves are significant) 
cannot be calculated, if we know the expression for only. To calculate this field one 
or another version of the perturbation theory must be used. 

We suggest to use Eq.(29) as a verification that the applied perturbation theory is 
accurate. 



5 CONCLUSIONS 




(29) 
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In some sense the calculation of the ohmic losses in the p-polarization state is much 
more simpl. The point is that if inequalities (1) are fulfilled the magnetic field penetrates 
into the grooves on the surface. In this case in the main approximation the magnetic field 
in the vicinity of the grooves is the same as near a flat surface of the perfect conductor. 
Thus the entire inner surface of the grooves absorbs the electromagnetic wave. In M for 
an arbitrary one-dimensional periodic surface this argumentation allowed us to obtain 
rather simple expression for the leading terms of the effective surface impedance in 
the p- polarization state: 

g =ReC^--i^f. (30) 
Se/ s 2b 2b K J 

In Eq.(30) 2b is the period of the surface structure, L is the length of the surface profile 

per one period and S is the area of the groove in the plane x\X^. If the effective impedance 

is known, there are no problems with calculation of the ohmic losses: 

Qf = ^\H \ 2 ReQ. (31) 

However, some difficulties can arise when the perturbation theory is applied to calculate 
the evanescent waves near rough surface in the p- polarization state. In the present 
publication we do not examine this question. 
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6 APPENDIX 

Let us demonstrate that the solution X^^kb) (I = 1,2 :) of the matrix Eq.(18) guaran- 
tees the fulfillment of the energy conservation law for an arbitrary value of kb < 1. The 
equation (18) corresponds to the infinitely conducting grating (£ = 0), and, consequently, 
the reflection coefficient has to be equal to one. With regard to our notations this means 
that substituting X^fL^kb) in Eq.(15) we have to verify that 1 — R = 1 — |e^~ f = 0. When 
C = 0, the explicit form of 1 — R is 

*h (2/-1) 2 ' W l ^L 1 {2l-lf{2m-lf 

(compare with Eq.(24); /i = 2a/ b). When writing Eq.(A.l) we take into account that 
W41 = 1/(2/ - l) 2 (see Eq.(14.2)). 

We seek X^l-iikb) in the form of a power series 

oo 
9=0 

Then it is convenient to rewrite Eq.(18.1) in the form 

oo i [ oo ] 
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The elements of the matrix $2r-i,2m-i depend on (kb) 2 only, 



$2r-i,2m-i = (2m - 1) £ <W 1 ~ (kb/nqf £>2-i,2m-i + ^ (0 ^ 1 - ^(2^- 1) ) 6rm > 

(A.Sb) 

where D$_ 12m _ x is defined by Eq.(18.3). 

Let us present $2r-i,2m-i as a series expansion: 

oo 

$ 2 ,-i, 2m -i = E *2Y P -Wi(M>) 2p , (A4) 

p=0 

When we substitute Eq.(A.2) in Eq.(A.3a) and set equal the terms, corresponding to 
the same powers of (kb), we obtain a set of equations relating the coefficients l^m-i 
with different values of the superscripts. Let us examine the obtained equations in the 
consecutive order. For q = (the term independent of kb) we have 

Evidently, Y^ ^ = A^L^O). Note, the coefficients Y 2 ^- X being solutions of the linear set 
of equations with real parameters, are real numbers. 
Next, for q > 1 we have 

Sq °° (2 ) ( 2 ) 2 /i 00 Y2 ^ x 

The first sum in the left-hand side of Eq.(A.5b) is taken over p < s q , where s q — s — 1, if 
q — 2s — 1, and s g = s, if g = 2s. 

The structure of the right-hand sides of Eqs.(A.5b) suggests to seek the coefficients 

*2?-i (? > 1) in the form 

Ylt, = YSl^ + Z^. (A6) 

Then for g = 1 we have 

4^1=0; Y 2 % = y£1 x T . (A.7a) 

We see that the coefficients Y 2 ^_ x are imaginary numbers. At the next step (q = 2) the 
equation for Z^J_ X is 



X/ ^2r-l,2m-1^2r-l + X/ ^2r-l,2m-1^2r-l _ 0- (A. 76) 



(0) 7 (2) v-^) v (0) 

^2r-l,2m-l Z/ 2r-l T 
r=l r=l 

With regard to Eq.(A.6) and Eqs.(A.7) it is evident that both Z 2r _ x and Y 2r _ x are real 
numbers. When continuing this procedure, we find out that all the coefficients Y^r-i^ 
are imaginary numbers and all the coefficients Y^-i are real numbers. On the contrary, 
T q is a real number when q is odd, and T q is an imaginary number when q is even. 

In the consecutive order we wrote down equations for Z 2 f_ x . Examining these equa- 
tions we obtained formulae relating Z^J_ X with the other coefficients z!fj_ x with r < q. At 
the last step we eliminated the functions Z 2 }_ x with the aid of Eq.(A.6). 
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Here, omitting the intermediate calculations, we present only the recurrence formula 
for the coefficients Y^r-i ^ obtained with the aid of the aforementioned procedure: 

ifir" = &t t^-t? e" i^.d«r I> -")'. (As) 

27r s =i (2s - 1) p= o 

To obtain Eq.(A.8) we used the definition of the sum T q (see Eq.(A.5b)). We do not need 
the explicit form of Y^r-i to prove that the right-hand side of Eq.(A.f) vanish. 

Now let us write the expression for 1 — R in the form of a series in powers of kb with 
( r > 1) being the coefficients of this series. It is easy to see that since the coefficients 
Y^r-i are re& l numbers, v-y — 0. If r > 1, we have 

,, oo v( r_1 ) ('v( r_1 )'l* /,,\2 oo 1 r-2 

vr^ (2/ — 1) U; ^ 1 (2/-l) 2 (2m-l) 2 ^ 2 ™ 1 l2 ^ 

(A9) 

Let us show that v r = for an arbitrary r. Indeed, if r = 2s — 1, taking into account 
that all the coefficients Y^* -1 ^ are real numbers, we see that the first term in the right- 
hand side of Eq.(A.9) is equal to zero. In the second term decomposing the sum over q 
into two sums over odd and even values of q, after a simple transformation we obtain 



^.-i = -(z-) 2 E T^r-^ — ^eV^t 9) - 3) + (Y£r*-*)\ (aio) 

W i, m =i (2Z — 1) (2m- 1) q=0 



and, consequently, ^2s-i = 0. 

Next, let us examine the case r = 2s. If we take into account that all the coefficients 
are imaginary numbers, from Eq.(A.9) we obtain 



,, oo V (2r-1) 2 oo 1 2(r-l) 

tt^(2/-I) 2 U; ^ (2Z - l) 2 (2m - l) 2 ^ 2m 1 2 '- lJ 



Now we make use of Eq.(A.9). From this equation it follows that 



(Alia) 



oo v (2s-l) . oo 2(s-l) 

fe(2/-l) 2 2vr^ 1 (2/-l) 2 (2m-l) 2 ^ 2 ™ U ^ J 1 J 

From the last equation we immediately have u 2s = 0. 

Thus we have shown that if ( = and kb < 1, the solution of the matrix equation 
(18.1) provides the fulfillment of the equality 1 — R = in agreement with the energy 
conservation law. 
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LIST OF FIGURES 

Fig.l. The grating configuration; 2b is the period, 2a is the width of the grooves 
throats. 

Fig. 2. The dimensionless values q s = ttQ^ ut ^ /cb\E \ 2 (solid line) andgj = ttQ^ / cb\E \ 2 
(dashed line) as functions of the ratio a/6. The crosses show the values of q v = nQ^ / cb\E \ 2 
calculated with the aid of Eq.(26). 

Fig.3. The ratios Nr = Re</?2n-i(0)/v4n-i(0) (solid lines) and Ni = Im<^2n-i(0)/v2n-i(0) 
(dashed line) in the point x\ — for two values of (5: (3 = 0.05 (Fig.3a) and (3 = 0.01 
(Fig.3b). 
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